Multiplier phenomenology in random multiplicative cascade processes 
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We demonstrate that the correlations observed in con- 
ditioned multiplier distributions of the energy dissipation in 
fully developed turbulence can be understood as an unavoid- 
able artefact of the observation procedure. Taking the latter 
into account, all reported properties of both unconditioned 
and conditioned multiplier distributions can be reproduced 
by cascade models with uncorrelated random weights if their 
bivariate splitting function is non-energy conserving. For the 
a-model we show that the simulated multiplier distributions 
converge to a limiting form, which is very close to the experi- 
mentally observed one. If random translations of the observa- 
tion window are accounted for, also the subtle effects found in 
conditioned multiplier distributions are precisely reproduced. 
PACS: 47.27.Eq, 05.40. +j, 02.50.Sk 

Random multiplicative cascade processes frequently 
serve as phenomenological models for the study of a vari- 
ety of complex systems exhibiting multifractal behaviour. 
In particular the energy dissipation field of fully devel- 
oped turbulent flows has been most successfully modelled 
in such terms, physically motivated by Richardson's pic- 
ture of energy transfer from large to small scales by ran- 
dom breakups of eddies [Q . 

Binary multiplicative cascade models relate the energy 
flux El at some integral scale L to Er = E^Wi ■ ■ ■ Wj 
contained in a subinterval of size r = 2~'^ at scale J 
by a product of mutually independent random weights 
Wj. More precisely, the energy flux Ej^ , contained in 
the interval k with length r = L/2-' splits into a left 
(L) and right (R) offspring interval, each of length r/2, 
whereby the content propagates according to E^^'^'^ — 
WlE\^^ and = WrE^s^\ respectively. For each 

breakup the two random weights Wl.r > are chosen, 
independently from any preceding breakup, according to 
a joint probability density p(Wl, Wr) with expectations 
{Wl,r)=1/2. The latter we denote as splitting function. 
Note that for a full description of a binary breakup the 
specification of the joint density is necessary. For the 
special case where the splitting function is concentrated 
along the diagonal Wl + Wr = 1, i.e. p{Wl,Wr) = 
p{Wl)S{Wl + Wr — 1), each breakup strictly conserves 
energy; for all other forms, however, the relation E^i^^^ + 

E^-'^^\' = El" holds only on average. 

A simple example is the splitting function 

p{Wl,Wr)^ (1) 
[piSiWL - (1 - a)/2)+p25{WL - (1 +/3)/2)] * 



[piS{Wr -{1- a)/2) +P25{Wr - (1 +/5)/2)] , 

which is known as the a-model Q . The probabilities pi = 
l—p2 — (3/(a+(3) are determined by energy conservation 
in the mean: pi{l — a)+p2{l+P) = 1. In contrast to the 
otherwise similar p-model suggested in Q, the a-model 
does not conserve energy in each local splitting, since 
with probability p^ + p\ we have Wl + Wr, ^ 1. 

Multiplicative cascade models are based on two as- 
sumptions: (a) the existence of a scale-independent split- 
ting function p{Wl,Wr) and (b) statistical indepen- 
dence of the random weights Wl,r. at one breakup from 
those of any other breakup. Once p{Wl , Wr) is chosen or 
deduced from experiment, the assumptions (a) and (b) 
allow to determine all moments and scaling exponents 
of the energy dissipation Sr—Er/r by inverse Laplace 
transforms Many experimental measurements of 

moments and scaling exponents confirm that this sim- 
ple construction reproduces the measured multifractal as- 
pects of the energy dissipation field amazingly well (e.g. 




A series of experimental investigations aim at a 
direct study of the random weights W by analysing the 
distributions of closely related observables, so-called mul- 
tipliers M (or breakup coefficients q—2M in the case of 
binary breakups), whose operational definition is given 
below. These studies reveal that assumption (a) is true 
over a decade of scales, with the restriction that the scale- 
invariant distribution p{M) depends also on the relative 
position of the offspring to the parent interval; however, 
assumption (b) is clearly violated. Significant correla- 
tions between multipliers at adjacent scales are observed. 
Such correlations apparently obscure the validity of un- 
correlated multiplicative cascade models for fully devel- 
oped turbulence, and more elaborate models, such as the 
'correlated p-model' [||, were suggested. These experi- 
mental findings gave already rise to a critical discussion 
of the limitations of multiplier phenomenology [p^ . In 
this Letter we present a further clarification of the exper- 
imental results. 

In order to explain the observed multiplier correla- 
tions within the framework of binary multiplicative cas- 
cade processes two considerations are made: (i) the bi- 
variate splitting function p{Wl,Wr) should not be as- 
sumed to be energy conserving; this is justified by the 
experimental restriction to measure the energy dissipa- 
tion field of three-dimensional turbulence only along a 
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one-dimensional cut (i.e. from the velocity time series 
obtained from anemometers and employment of Taylor's 
frozen flow hypotheses). Consideration (ii) concerns the 
obvious non-homogeneity of cascade processes. Due to 
their hierarchical structure the n-point correlation func- 
tions of the generated energy dissipation field arc not 
translationally invariant (cf . pH | for a visualisation) ; mul- 
tipliers in dyadic intervals will be different if the obser- 
vation window L is shifted by an arbitrary amount, say 
Aa;. In real world experiments, however, the observation 
window is placed in no relation with the unobservable 
hierarchical structure of the cascade, so that implicitly 
an averaging over uniform random translations is per- 
formed. The effect of such random translations was hith- 
erto assumed to be negligible. We test this assumption by 
introduction of random shifts before calculating the mul- 
tipliers. It is especially this latter operation that leads to 
a full explanation of all available experimental findings of 
conditional multiplier distributions. In the following we 
discuss the implications of above considerations (i) and 
(ii) one after the other. 

It has already been pointed out in ||l^ that the non- 
conservation of energy in the splitting function leads to 
deviations from perfect multifractal scaling. Here we fo- 
cus on its influence on multiplier distributions obtained 
from simulated realizations e'^^ after J cascade steps. 
Depending on the relative position of parent to offspring 
intervals, left, right and centred multipliers at scale j and 
position k are operationally defined as 
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and the multiplier distributions p{M^^^) are obtained by 
histogramming all < fc < 2-' possible multipliers at a 
given scale. 

It is important to note that the experimentally mea- 
sured 'backward' energies 
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are obtained by successive summations from the finest 
resolution scale J to larger ones and are generally not 
equal to the 'forward' energy densities Ej: , which arise 
as intermediate states in the evolution of the cascade 
from larger to smaller scales. For a clearer distinction 
we denote the former with a bar. Since, by definition, 
the multipliers satisfy M^.l + Mk^n. = 1, they effectively 
enforce local energy conservation and therefore give only 
incomplete information on the true parent-offspring rela- 
tion of non-energy conserving breakups. 



For the special case of energy-conserving splitting func- 
tions, such as the p-model [||, the multipliers (||) do give 
a faithful representation of local breakups. Here, the left 
and right multiplier distributions are equal to the split- 
ting function p{W) along the Wl + Wr = 1 diagonal and 
do not depend on the scale j. 

The situation is quite different for non-energy conserv- 
ing models, such as the a-model. In this case there is 
no simple analytical relation between the splitting func- 
tion (|l|) and the resulting distributions of multipliers (0); 
in fact, the mapping from (|l|) to p{M^-^^) via Eqs. ^) 
and (^ acts effectively as a non-linear smoothing oper- 
ation. This we demonstrate in the following simulation 
results: The scale dependent left|^ multiplier distribution 
p(Af|-'-*) is depicted in Fig. 1, where 10^ configurations 
were generated, each with J = 9 cascade steps. While 
the weights Wl,r in the forward evolution may take only 
two values (1 — a)/2 and (1 -|- /9)/2, the multipliers (|^) 
take more and more distinct values as the scale differ- 
ence J — j increases and becomes soon quasi-continuous. 
Already after three backward steps the left (and right) 
multiplier distributions apparently converge to a limit- 
ing form which comes surprisingly close to a parameter- 
isation of the experimentally observed multiplier distri- 
bution shown as continuous line. The parameters 
used in this simulation were a = 0.3, (3 = 0.65, but this 
finding also holds for other parameters as, for example, 
a — j3 — 0.4 or a = 0.5, /3 = 0.3; all these choices have 
approximately the same second order splitting moment 
{WD = /o p{Wl, Wb) Wl dWLdWR. 

A similar convergence follows for centred multiplier 
distributions p{Mc) with the distinction that the corre- 
sponding limiting density appears to be narrower (shown 
in Fig. 2a for j = 3) compared to p{Ml); this latter fea- 
ture is quite in agreement with the observations in 
Analogous results are found for multipliers generalised 
to arbitrary length ratios between offspring and parent 
intervals. In Fig. 2a the left-skewed histograms show cor- 
responding limiting distributions for left and central mul- 
tipliers for two cascade steps j ~> j + 2 and thus relate 
intervals with length ratio 1/4. 

In above comparisons the unavoidable random transla- 
tions of the observation window with respect to the cas- 
cade position (as discussed in suggestion (ii) above) has 
not been accounted for. Following the procedure intro- 
duced in Ref . , we adopt a scheme of random transla- 
tions in the following way: for the target resolution scale 
J of a given integral length scale L a longer cascade re- 
alization with J + 3 steps is generated (corresponding to 



^In all simulations the distributions for left (L) and right 
(R) multipliers are identical; therefore only one (L) is shown. 
Differences occur, however, with centred multipliers and the 
latter are shown in Figs. 2 and 3. 
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an integral length scale 8L) and an observation window 
of size L is shifted randomly by t bins within 8L. In 
other words, only the 2"^ bins of the generated eI.'!~^^\ 
which lie within the randomly placed observation win- 
dow are considered, giving a simulated and translated 
Elf^ = E\j^f\ where t is a uniformly distributed inte- 
ger within [0, 7-2'^ — 1]. Then the multipliers are deter- 
mined again by and sampled over a large number of 
configurations and random shifts t. 

There is no dramatic change in unconditioned multi- 
plier distributions after addition of random translations; 
the small effect is illustrated in Fig. 2b compared to 
Fig. 2a. The left multiplier distribution p{Ml) now is 
even closer to the experimentally deduced scale-invariant 
parametrisation of Ref. 1^; this property also holds for 
the other two parametrisations used, {a, (3) = (0.5, 0.3) 
and (0.4,0.4). Note, that the centred multiplier distribu- 
tion p{Mc) now suffers a small but noticeable asymme- 
try which is also seen in experimental observations 
Moreover, the gaps close to the endpoints at M = and 
M = 1 are now filled, which is a point of interest in 
the discussion of Novikov's 'gap-theorem' and 
weakens any conclusions drawn from it^. 

So far the multiplier distributions we have looked at 
are unconditioned. In the experimental analyses 
conditioned multiplier distributions of the form 
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have been shown; they correlate a parent multiplier 
with the multiplier Af^"* of its offspring, where 
A stands for L,C,R respectively. 

The addition of random translations introduces signif- 
icant correlations among multipliers at different scales, 
which are well reflected in the conditioned distributions 
(^. Figs. 3a and 3b illustrate the correlations between 
offspring and parent multiplier for the centred and left 
case respectively. Compared to the unconditioned den- 
sity p{M^J'') the conditioned one p(M^^^|0 < M^^"^^ < 
1/2) is skewed to the left and p(Ar^^^|l/2 < M^J~'^^ < 1) 
is skewed to the right; the centred offspring and parent 
multipliers are positively correlated. Moreover, the max- 
imum value of the former is higher than the later. This 
result is identical to the experimental finding of Ref. . 

For the conditioned left multiplier distributions 
p(Mj;'^|Mmin < M^^^^^ < Mmax) we show results in 
Fig. 3b; again the parameters {a,f3) — (0.3,0.65) for 
the ck-model have been used. The conditioning on the 
sub-range [Afmin , Afmax] = [0.2,0.4] leads to narrowing 



of the distribution and conditioning on [Afmin, Afmax] = 
[0.6, 0.8] to a broadening. Again, this is in perfect agree- 
ment with the experimental finding in Ref. [^ . 

We report without figures that for the parameter 
choice {a, f3) = (0.4,0.4) the effects seen in Fig. 3b 
become weaker and vanish almost completely for the 
choice (a,/3) = (0.5,0.3). Also Fig. 3a is modified: 
for {a, (3) — (0.5,0.3) the maximum of the left skewed 
distribution (1) is lower than the one of (2), while for 
(a,/J) = (0.4,0.4) they are about equal. A noticeable 
left /right shift, however, remains in all three cases. 

From these observations we come to the following con- 
clusion: with non-energy conserving splitting functions 
and the inclusion of random translations the uncondi- 
tioned multiplier distributions observed in the data [||j|] 
are reproduced quite naturally; moreover, once the in- 
put splitting function is skewed in a certain direction 
also the correct conditioned multiplier distributions are 
deduced. Very similar numerical observations are also 
obtained for different choices of splitting functions [Q. 
Since the simulated cascades satisfy both assumptions 

(a) and (b) above, we regard the comparable violation of 

(b) in experiments and Fig. 3 as an unavoidable artefact 
of the observation procedure. 

Related to the multiplier phenomenology is the prob- 
lem of how to extract the correct (multifractal) scaling 
exponents. The findings of Ref. favour the left /right 
over the central multipliers, but there the random trans- 
lations were not accounted for. Certainly the inclusion 
of the latter modifies the scaling exponents to some ex- 
tent [|l^] , the more noticeable the higher the order of the 
underlying moments. Due to these operational ambigu- 
ities we feel that it should not be a matter of how to 
extract scaling exponents, but how to deduce the opti- 
mal non-energy conserving and skewed splitting function 
from data. A starting point are the findings of Ref. [^ , 
which show how to extract bivariate splitting functions; 
however, the aspect of translation invariance has not been 
taken into account. A satisfactory solution of this intri- 
cate inverse problem will only be feasible once additional 
observables are studied experimentally, such as n-point 
correlation functions or their wavelet compressed form 
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^The input splitting function of the a-model has gaps (1 — 
a)/2 at W=0 and (1 — /3)/2 at W—l; both gaps disappear in 
a multiplier analysis. 
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FIG. 1. Convergence of the left multiplier distributions p{Mf^ ) of the non-energy conserving a-model (|l|) to a 
quasi-continuous limiting form. Parameters are J = 9, a = 0.3 and /3 = 0.65. For comparison the experimentally de- 
duced Beta-function parametrisation of the left (right) unconditioned multiplier distribution with /3 = 3.2 ^ is shown as a 
solid line. 
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FIG. 2. (a) Left (dashed) and centred (solid) multiplier distributions for the a-model (|l|). For the one step multipliers 
(j ~* J + 1) the length ratio between offspring and parent intervals is A = i while for the two-step multipliers (j ^ j + 2) 
we have \ — j. Parameters are J = 9, q = 0.3 and (3 = 0.65.(b) Same as (a), but with random translations applied 
before calculating the multipliers. For comparison the experimentally deduced Beta-function parametrisation of the left (right) 
multiplier distribution (j j + 1) with /3 — 3.2 M is shown as a thick solid line. 
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FIG. 3. Conditioned multiplier distributions for the a-model, satisfying strictly conditions (a) and (b) in the text. Parameters 
are J — 9, a — 0.3, (3 = 0.65 and random translations were applied before calculating the multipliers, (a) Central distributions 
conditioned on central parent multipliers: (1) p(M^^'|0 < M^^"^' < 1/2) and (2) p{M^J''\l/2 < M^J''^'' < 1). (b) Left 
distributions conditioned on left parent multipliers: (1) p{M[^^ \0.2 < M^f"^' < 0.4) and (2) p(M}f'' |0.6 < M^-'"^' < 0.8) 
compared to the Beta-function parametrisation of unconditioned distributions in Fig. 1 and Fig. 2 (full line). 
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